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In attempting to understand the endogenic processes which have shaped the surface of an icy satellite, it 
is desirable to quantify the failure strength of the satellite’s lithosphere. In a crust that is fractured on a 
large scale, frictional sliding along pre-existing fractures occurs in response to lower differential stresses 
than required to initiate fracture of pristine rock, thus governing failure of a brittle lithosphere. Failure is 
predicted along favorably -oriented fracture planes; if fractures of all orientations are assumed to be present in 
the crust (as is expected of a heavily cratered lithosphere), frictional failure relations are directly applicable. 
The Coulomb criterion predicts that the shear stress a T and normal stress a n components on a fracture plane 
at failure are related as o t = pa n + S Q where S 0 is the cohesion and p is the coefficient of friction. At 
moderate to high pressures, the frictional sliding strength of most materials is found to be a x = 0.85<r rt 
[1]. This strength is largely independent of rock type and strength, strain rate, and temperature. 

The frictional failure curve for ice, however, deviates from this law [2]. At low temperatures (77 < T < 
115 K) and moderate pressures ( a n < 200 bar), ice frictional failure can be described by cx T = 0.55cr rt + 10 
bar, independent of temperature or sliding rate. The frictional failure curve for real surfaces probably does 
not have a cohesion intercept, but may instead have vertical tangency to the shear stress axis, conceptually 
intercepting the origin [3]. Constraining the frictional failure line for ice to pass through the origin, cr T = 
0.69cx n is obtained [2]. 

At low confining pressures (oj 5 30 bar), frictional sliding of a variety of rock types is strongly 
dependent on the surface roughness of the pre-existing fractures [1,3]. Such low stresses are applicable in 
considering frictional failure of low gravity satellites such as Miranda. The frictional strength of ice has 
been tested down to oj = 3 bar [2], but the effect of surface roughness has not been directly evaluated. 

Near the horizontal surface of a planet, lithospheric stresses may be resolved into three principal stress 
directions, one of which is vertical and the other two are horizontal. In a lithosphere of average density p 
subject to gravitational acceleration g , the vertical principal stress can be equated (in the absence of pore 
pressure) to the lithostatic pressure pgz at some depth z. In the case of horizontal compression, the 
maximum compressional stress <jj is horizontal and the minimum compressional (maximum tensional) 
stress Oj is vertical; in the case of horizontal extension, the opposite is true. With one principal stress thus 
known, the other necessary to induce failure may be determined by 
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for some value of friction and negligible cohesion [4]. 

The maximum horizontal stress that may be supported in a lithosphere may be represented as the 
difference o H - c v between the horizontal and vertical principal stresses [5], with compression being 
considered positive. Frictional failure stress is plotted a function of depth in Figure 1 for a crustal density 
of 1000 kg m'3, gravity = 0.09 m s" 2 , and friction values n = 0.5 and 0.7, chosen as bounds to the ice 
frictional failure data. 

As illustrated by Figure 1, frictional failure strength of the brittle lithosphere increases with depth as 
lithostatic pressure increases for both compression (crj vertical, positive stress) and extension (ct; vertical, 
negative stress) until the brittle-ductile transition depth is reached. Below this depth, lithospheric strength 
is controlled by creep of the ductile crust and decreases with depth. Behavior of the ductile lithosphere is 
dependent upon composition, temperature, and strain rate. Plotted on Figure 1 are curves that represent the 
ductile strength of cold ice based on the data of [6], A strain rate of 10‘ 16 s' 1 and a surface temperature of 
60 K are chosen, and the labeled curves represent thermal gradients of 2, 5, 10, and 20 K km' 1 . 
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The plot shows that the maximum strength of Miranda's lithosphere (that at the brittle-ductile 
transition) is small. Failure in compression (reverse faulting) is induced by < 50 bar of horizontal stress, a 
stress close to that which may cause folding of an ice crust [7], indicating that both deformation styles 
might be active in a compressional regime. Less than 20 bar could cause failure in extension (normal 
faulting), a stress which could be supplied by satellite expansion of < 0.05% [8]. 

If troughs in Miranda's Elsinore Corona are graben whose faults intersected near the brittle-ductile depth 
transition at the time of their formation, a transition depth of ~5 to 10 km is implied [9]. This necessitates 
a paleo-thermal gradient -10 to 20 K km 1 . A steeper thermal gradient (and/or a lesser strain rate) is 
suggested if narrow troughs (-1 km wide) in Elsinore and Inverness Coronae are graben that formed in a 
similar fashion. A lower thermal gradient earlier in Miranda’s history may have facilitated formation of 
broader troughs, which appear to be more ancient [10]. 

Comparing the results of Figure 1 to those for Ganymede [2], it is apparent that lower satellite gravity 
results in much lower frictional strength and in a somewhat deeper brittle lithosphere. Lower surface 
temperature significantly lowers the depth to ductile ice behavior for a given lithospheric stress, thermal 
gradient, and strain rate. The results here can be directly applied to Enceladus, which has a gravity and a 
surface temperature similar to those of Miranda. Enceladus displays subparallel troughs ~4 km wide that 
may be graben which formed in a brittle lithosphere -3 km deep. For a water ice lithosphere, low stress 
but high paleo-thermal gradient (>20 K km" 1 ) is indicated. Such a high gradient is likewise necessary to 
account for the shapes of some craters on Enceladus if a water ice lithosphere is assumed [1 1]. 

The results of Figure 1 are somewhat suspect because of poor constraints on the creep of ice at very 
low temperatures. In addition, the lithospheres of Miranda and Enceladus are probably not composed of 
pure water ice: Miranda possesses a dark contaminant, and both lithospheres may be composed of a mixture 
of water and other ices, which could allow ductile behavior at lower temperatures and more gentle thermal 
gradients. In addition, the possible effect of fracture roughness needs to be evaluated, as this might cause 
the friction coefficient p to vary be more than an order of magnitude, greatly affecting the slope of the 
frictional strength line. 

The low frictional strength of low-gravity satellites permits relatively minor lithospheric stress to 
result in surface defamation; therefore, only small horizontal stresses are necessary to account for observed 
endogenic deformation on Miranda and Enceladus. The low surface temperature of these satellites 
necessitates steep thermal gradients to produce shallow ductile behavior of a water ice lithosphere. 



Figure 1. Maximum 
horizontal stress which can 
be supported in an ice 
lithosphere for Miranda (or 
Enceladus). Strength of the 
brittle lithosphere is 
controlled by frictional 
failure and is bracketed by 
straight lines for n = 0.5 and 
0.7 for gravity g * 0.09 
m s" 2 . Ductile strength 
curves for cold ice are for a 
strain rate of 10" 16 s* 1 , a 
surface temperature of 60 K, 
and thermal gradients as 
indicated. Dots indicate 
depths at which T = 195 K is 
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attained, below which the 
low-temperature ice flow law 
parameters used here are no 
longer strictly valid [6]. 
Compressional stress is 
taken to be positive. 
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